home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_hpmv.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  11.5 KB  |  380 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_hpmv (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int order = 101;
  13.    int uplo = 121;
  14.    float alpha[2] = {-0.3f, 0.1f};
  15.    float beta[2] = {0.0f, 1.0f};
  16.    int N = 2;
  17.    float A[] = { 0.339f, -0.102f, 0.908f, 0.097f, -0.808f, 0.236f };
  18.    float X[] = { 0.993f, -0.502f, -0.653f, 0.796f };
  19.    int incX = -1;
  20.    float Y[] = { -0.35f, 0.339f, -0.269f, -0.122f };
  21.    int incY = -1;
  22.    float y_expected[] = { -0.0627557f, -0.839323f, -0.0877262f, -0.169208f };
  23.    cblas_chpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  24.    {
  25.      int i;
  26.      for (i = 0; i < 2; i++) {
  27.        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chpmv(case 1118) real");
  28.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chpmv(case 1118) imag");
  29.      };
  30.    };
  31.   };
  32.  
  33.  
  34.   {
  35.    int order = 101;
  36.    int uplo = 121;
  37.    float alpha[2] = {-0.3f, 0.1f};
  38.    float beta[2] = {0.0f, 1.0f};
  39.    int N = 2;
  40.    float A[] = { 0.339f, -0.102f, 0.908f, 0.097f, -0.808f, 0.236f };
  41.    float X[] = { 0.993f, -0.502f, -0.653f, 0.796f };
  42.    int incX = -1;
  43.    float Y[] = { -0.35f, 0.339f, -0.269f, -0.122f };
  44.    int incY = -1;
  45.    float y_expected[] = { -0.0627557f, -0.839323f, -0.0877262f, -0.169208f };
  46.    cblas_chpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  47.    {
  48.      int i;
  49.      for (i = 0; i < 2; i++) {
  50.        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chpmv(case 1119) real");
  51.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chpmv(case 1119) imag");
  52.      };
  53.    };
  54.   };
  55.  
  56.  
  57.   {
  58.    int order = 101;
  59.    int uplo = 122;
  60.    float alpha[2] = {-0.3f, 0.1f};
  61.    float beta[2] = {0.0f, 1.0f};
  62.    int N = 2;
  63.    float A[] = { 0.339f, -0.102f, 0.908f, 0.097f, -0.808f, 0.236f };
  64.    float X[] = { 0.993f, -0.502f, -0.653f, 0.796f };
  65.    int incX = -1;
  66.    float Y[] = { -0.35f, 0.339f, -0.269f, -0.122f };
  67.    int incY = -1;
  68.    float y_expected[] = { -0.0037603f, -0.816761f, -0.0392456f, -0.121154f };
  69.    cblas_chpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  70.    {
  71.      int i;
  72.      for (i = 0; i < 2; i++) {
  73.        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chpmv(case 1120) real");
  74.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chpmv(case 1120) imag");
  75.      };
  76.    };
  77.   };
  78.  
  79.  
  80.   {
  81.    int order = 101;
  82.    int uplo = 122;
  83.    float alpha[2] = {-0.3f, 0.1f};
  84.    float beta[2] = {0.0f, 1.0f};
  85.    int N = 2;
  86.    float A[] = { 0.339f, -0.102f, 0.908f, 0.097f, -0.808f, 0.236f };
  87.    float X[] = { 0.993f, -0.502f, -0.653f, 0.796f };
  88.    int incX = -1;
  89.    float Y[] = { -0.35f, 0.339f, -0.269f, -0.122f };
  90.    int incY = -1;
  91.    float y_expected[] = { -0.0037603f, -0.816761f, -0.0392456f, -0.121154f };
  92.    cblas_chpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  93.    {
  94.      int i;
  95.      for (i = 0; i < 2; i++) {
  96.        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chpmv(case 1121) real");
  97.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chpmv(case 1121) imag");
  98.      };
  99.    };
  100.   };
  101.  
  102.  
  103.   {
  104.    int order = 102;
  105.    int uplo = 121;
  106.    float alpha[2] = {-0.3f, 0.1f};
  107.    float beta[2] = {0.0f, 1.0f};
  108.    int N = 2;
  109.    float A[] = { 0.339f, -0.102f, 0.908f, 0.097f, -0.808f, 0.236f };
  110.    float X[] = { 0.993f, -0.502f, -0.653f, 0.796f };
  111.    int incX = -1;
  112.    float Y[] = { -0.35f, 0.339f, -0.269f, -0.122f };
  113.    int incY = -1;
  114.    float y_expected[] = { -0.0627557f, -0.839323f, -0.0877262f, -0.169208f };
  115.    cblas_chpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  116.    {
  117.      int i;
  118.      for (i = 0; i < 2; i++) {
  119.        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chpmv(case 1122) real");
  120.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chpmv(case 1122) imag");
  121.      };
  122.    };
  123.   };
  124.  
  125.  
  126.   {
  127.    int order = 102;
  128.    int uplo = 121;
  129.    float alpha[2] = {-0.3f, 0.1f};
  130.    float beta[2] = {0.0f, 1.0f};
  131.    int N = 2;
  132.    float A[] = { 0.339f, -0.102f, 0.908f, 0.097f, -0.808f, 0.236f };
  133.    float X[] = { 0.993f, -0.502f, -0.653f, 0.796f };
  134.    int incX = -1;
  135.    float Y[] = { -0.35f, 0.339f, -0.269f, -0.122f };
  136.    int incY = -1;
  137.    float y_expected[] = { -0.0627557f, -0.839323f, -0.0877262f, -0.169208f };
  138.    cblas_chpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  139.    {
  140.      int i;
  141.      for (i = 0; i < 2; i++) {
  142.        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chpmv(case 1123) real");
  143.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chpmv(case 1123) imag");
  144.      };
  145.    };
  146.   };
  147.  
  148.  
  149.   {
  150.    int order = 102;
  151.    int uplo = 122;
  152.    float alpha[2] = {-0.3f, 0.1f};
  153.    float beta[2] = {0.0f, 1.0f};
  154.    int N = 2;
  155.    float A[] = { 0.339f, -0.102f, 0.908f, 0.097f, -0.808f, 0.236f };
  156.    float X[] = { 0.993f, -0.502f, -0.653f, 0.796f };
  157.    int incX = -1;
  158.    float Y[] = { -0.35f, 0.339f, -0.269f, -0.122f };
  159.    int incY = -1;
  160.    float y_expected[] = { -0.0037603f, -0.816761f, -0.0392456f, -0.121154f };
  161.    cblas_chpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  162.    {
  163.      int i;
  164.      for (i = 0; i < 2; i++) {
  165.        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chpmv(case 1124) real");
  166.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chpmv(case 1124) imag");
  167.      };
  168.    };
  169.   };
  170.  
  171.  
  172.   {
  173.    int order = 102;
  174.    int uplo = 122;
  175.    float alpha[2] = {-0.3f, 0.1f};
  176.    float beta[2] = {0.0f, 1.0f};
  177.    int N = 2;
  178.    float A[] = { 0.339f, -0.102f, 0.908f, 0.097f, -0.808f, 0.236f };
  179.    float X[] = { 0.993f, -0.502f, -0.653f, 0.796f };
  180.    int incX = -1;
  181.    float Y[] = { -0.35f, 0.339f, -0.269f, -0.122f };
  182.    int incY = -1;
  183.    float y_expected[] = { -0.0037603f, -0.816761f, -0.0392456f, -0.121154f };
  184.    cblas_chpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  185.    {
  186.      int i;
  187.      for (i = 0; i < 2; i++) {
  188.        gsl_test_rel(Y[2*i], y_expected[2*i], flteps, "chpmv(case 1125) real");
  189.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], flteps, "chpmv(case 1125) imag");
  190.      };
  191.    };
  192.   };
  193.  
  194.  
  195.   {
  196.    int order = 101;
  197.    int uplo = 121;
  198.    double alpha[2] = {1, 0};
  199.    double beta[2] = {0, 1};
  200.    int N = 2;
  201.    double A[] = { 0.543, -0.737, 0.281, -0.053, -0.098, -0.826 };
  202.    double X[] = { 0.67, -0.857, -0.613, -0.927 };
  203.    int incX = -1;
  204.    double Y[] = { -0.398, -0.934, -0.204, 0.183 };
  205.    int incY = -1;
  206.    double y_expected[] = { 0.745218, -0.60699, -0.37301, -0.983688 };
  207.    cblas_zhpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  208.    {
  209.      int i;
  210.      for (i = 0; i < 2; i++) {
  211.        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhpmv(case 1126) real");
  212.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhpmv(case 1126) imag");
  213.      };
  214.    };
  215.   };
  216.  
  217.  
  218.   {
  219.    int order = 101;
  220.    int uplo = 121;
  221.    double alpha[2] = {1, 0};
  222.    double beta[2] = {0, 1};
  223.    int N = 2;
  224.    double A[] = { 0.543, -0.737, 0.281, -0.053, -0.098, -0.826 };
  225.    double X[] = { 0.67, -0.857, -0.613, -0.927 };
  226.    int incX = -1;
  227.    double Y[] = { -0.398, -0.934, -0.204, 0.183 };
  228.    int incY = -1;
  229.    double y_expected[] = { 0.745218, -0.60699, -0.37301, -0.983688 };
  230.    cblas_zhpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  231.    {
  232.      int i;
  233.      for (i = 0; i < 2; i++) {
  234.        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhpmv(case 1127) real");
  235.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhpmv(case 1127) imag");
  236.      };
  237.    };
  238.   };
  239.  
  240.  
  241.   {
  242.    int order = 101;
  243.    int uplo = 122;
  244.    double alpha[2] = {1, 0};
  245.    double beta[2] = {0, 1};
  246.    int N = 2;
  247.    double A[] = { 0.543, -0.737, 0.281, -0.053, -0.098, -0.826 };
  248.    double X[] = { 0.67, -0.857, -0.613, -0.927 };
  249.    int incX = -1;
  250.    double Y[] = { -0.398, -0.934, -0.204, 0.183 };
  251.    int incY = -1;
  252.    double y_expected[] = { 0.646956, -0.542012, -0.282168, -0.912668 };
  253.    cblas_zhpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  254.    {
  255.      int i;
  256.      for (i = 0; i < 2; i++) {
  257.        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhpmv(case 1128) real");
  258.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhpmv(case 1128) imag");
  259.      };
  260.    };
  261.   };
  262.  
  263.  
  264.   {
  265.    int order = 101;
  266.    int uplo = 122;
  267.    double alpha[2] = {1, 0};
  268.    double beta[2] = {0, 1};
  269.    int N = 2;
  270.    double A[] = { 0.543, -0.737, 0.281, -0.053, -0.098, -0.826 };
  271.    double X[] = { 0.67, -0.857, -0.613, -0.927 };
  272.    int incX = -1;
  273.    double Y[] = { -0.398, -0.934, -0.204, 0.183 };
  274.    int incY = -1;
  275.    double y_expected[] = { 0.646956, -0.542012, -0.282168, -0.912668 };
  276.    cblas_zhpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  277.    {
  278.      int i;
  279.      for (i = 0; i < 2; i++) {
  280.        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhpmv(case 1129) real");
  281.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhpmv(case 1129) imag");
  282.      };
  283.    };
  284.   };
  285.  
  286.  
  287.   {
  288.    int order = 102;
  289.    int uplo = 121;
  290.    double alpha[2] = {1, 0};
  291.    double beta[2] = {0, 1};
  292.    int N = 2;
  293.    double A[] = { 0.543, -0.737, 0.281, -0.053, -0.098, -0.826 };
  294.    double X[] = { 0.67, -0.857, -0.613, -0.927 };
  295.    int incX = -1;
  296.    double Y[] = { -0.398, -0.934, -0.204, 0.183 };
  297.    int incY = -1;
  298.    double y_expected[] = { 0.745218, -0.60699, -0.37301, -0.983688 };
  299.    cblas_zhpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  300.    {
  301.      int i;
  302.      for (i = 0; i < 2; i++) {
  303.        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhpmv(case 1130) real");
  304.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhpmv(case 1130) imag");
  305.      };
  306.    };
  307.   };
  308.  
  309.  
  310.   {
  311.    int order = 102;
  312.    int uplo = 121;
  313.    double alpha[2] = {1, 0};
  314.    double beta[2] = {0, 1};
  315.    int N = 2;
  316.    double A[] = { 0.543, -0.737, 0.281, -0.053, -0.098, -0.826 };
  317.    double X[] = { 0.67, -0.857, -0.613, -0.927 };
  318.    int incX = -1;
  319.    double Y[] = { -0.398, -0.934, -0.204, 0.183 };
  320.    int incY = -1;
  321.    double y_expected[] = { 0.745218, -0.60699, -0.37301, -0.983688 };
  322.    cblas_zhpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  323.    {
  324.      int i;
  325.      for (i = 0; i < 2; i++) {
  326.        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhpmv(case 1131) real");
  327.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhpmv(case 1131) imag");
  328.      };
  329.    };
  330.   };
  331.  
  332.  
  333.   {
  334.    int order = 102;
  335.    int uplo = 122;
  336.    double alpha[2] = {1, 0};
  337.    double beta[2] = {0, 1};
  338.    int N = 2;
  339.    double A[] = { 0.543, -0.737, 0.281, -0.053, -0.098, -0.826 };
  340.    double X[] = { 0.67, -0.857, -0.613, -0.927 };
  341.    int incX = -1;
  342.    double Y[] = { -0.398, -0.934, -0.204, 0.183 };
  343.    int incY = -1;
  344.    double y_expected[] = { 0.646956, -0.542012, -0.282168, -0.912668 };
  345.    cblas_zhpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  346.    {
  347.      int i;
  348.      for (i = 0; i < 2; i++) {
  349.        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhpmv(case 1132) real");
  350.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhpmv(case 1132) imag");
  351.      };
  352.    };
  353.   };
  354.  
  355.  
  356.   {
  357.    int order = 102;
  358.    int uplo = 122;
  359.    double alpha[2] = {1, 0};
  360.    double beta[2] = {0, 1};
  361.    int N = 2;
  362.    double A[] = { 0.543, -0.737, 0.281, -0.053, -0.098, -0.826 };
  363.    double X[] = { 0.67, -0.857, -0.613, -0.927 };
  364.    int incX = -1;
  365.    double Y[] = { -0.398, -0.934, -0.204, 0.183 };
  366.    int incY = -1;
  367.    double y_expected[] = { 0.646956, -0.542012, -0.282168, -0.912668 };
  368.    cblas_zhpmv(order, uplo, N, alpha, A, X, incX, beta, Y, incY);
  369.    {
  370.      int i;
  371.      for (i = 0; i < 2; i++) {
  372.        gsl_test_rel(Y[2*i], y_expected[2*i], dbleps, "zhpmv(case 1133) real");
  373.        gsl_test_rel(Y[2*i+1], y_expected[2*i+1], dbleps, "zhpmv(case 1133) imag");
  374.      };
  375.    };
  376.   };
  377.  
  378.  
  379. }
  380.